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^ . We present an analysis of the decay D" K^n from FOCUS. From a sample of 234 

' events we find a branching ratio of pljoZ^jy-^+j ~ {0-'^30 ^ g'gg^ ± 0.031)% under the 

, assumption of no mixing and no CP violation. We also present limits on charm mixing. 

o 

^ , 

, 1. Introduction and event selection 

. Charm mixing remains elusive. In the Standard Model the charm mixing rate is 

I greatly suppressed by small CKM matrix elements and strong GIM suppression. 

■ This very fact, however, provides a unique opportunity to search for new physics. 

' , The search for charm mixing requires tagging the flavor of a neutral charm me- 

^ ' son at production and again tagging the flavor at decay. The production flavor is 

■ determined using Z)*+ D^T^t decays. The decay flavor is determined by recon- 
structing a Cabibbo favored (CF) decay with a charged kaon such as ^ K~ tt^ . 
A right-sign (RS) decay is defined as one in which the kaon and soft pion charges 
are opposite while same kaon and tTs charges are wrong-sign (WS) decays. The 

\ can also decay directly to X+tt" via a doubly Cabibbo suppressed decay (DCSD). 

The time dependent WS decay rate can be written as: 

Rws{t) = e-" {rd + VR^y'Tt + i (x'^ + y'') TH^^ (1) 

where the three terms correspond to DCSD, interference between mixing and DCSD, 
and mixing. Rd is the DCS branching ratio relative to the Cabibbo favored mode. 
The paramters x' = x cos 5k-k + y sin Sk-k and y' = y cos Sxn — x sin Sktt are rotated 
versions of the mixing parameters x = and y = ^ (mass and lifetime split- 
ting terms) and Sk-k is the strong phase between CF and DCS decay. Discovery of 
hadronic charm mixing requires separating the three components of Eq. fusing the 
different lifetime distributions. Note that Eq.^has no information on the sign of x' 
and therefore the relevent fit variable is x' . In this analysis, CP violation effects 
are not considered and charge conjugates are implied. 

FOCUS recorded data during the 1996-7 fixed-target run at Fermilab. A photon 
beam impinged on BeO targets. 16 silicon strip planes provide vertexing and track- 
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ing. Charged particles are tracked and momentum analyzed as they pass through 
up to two dipole magnets and up to five sets of multiwire proportional chambers. 

Three Cerenkov counters, two EM calorimeters, and two muon detectors identify 
particles. A candidate driven vertexing algorithm is used to reconstruct charm. For 
K~Tr~^, two tracks must verticize with CL > 2%. This candidate is used to 
seed the production vertex which must have CL > 1%. The soft pion 7r+ from the 
£)*+ — *• D^nf must be consistent with originating from the production vertex and 
the track is refit to the production vertex. Separating charm from hadronic back- 
ground is primarily accomplished by requiring the decay vertex be distinct from the 
production vertex. Other cuts are made to ensure the vertex is isolated and is more 
charm like than background like. Strong Cerenkov cuts are applied to the combina- 
tion of K~ and 7r+ especially if the invariant mass under a reflection of K~ — > tt" 
and ^K~^ is close to the D° mass. If more than one D*^ candidate is found for 
one D° candidate an additional cut may be made. If a RS (WS) candidate is found 
with 2.5 MeV < Q{D*) < 9.5 MeV then no WS (RS) candidates are allowed. 

2. Fit description and results 

A 3D binned maximum likelihood fit is used in this analysis. Two dimensions, 

M{D^) and Q{D*) separate signal from background while t{D^) distringuishes 
DCSD, mixing, and intcrcfcrcnce. The RS and WS data are fit simultaneously. 
The fit components generally have a shape determined by a Monte Carlo sim- 
ulation and a yield which is free to float, some with weak constraints imposed. 
The components making up the fit which arc modeled by the Monte Carlo arc 
RS signal, WS signal (DCSD, interference, and mixing), real D° K^n^ decay 
with a fake soft pion. reflections to both RS and WS {K~K^, 7r^7r+, Tr^Tr+Tr", 
ii'°7r-7r+), RS background {R-i+v and ii'-Tr+Tr"), and WS background (double 
misidentification of K^i^iy, K~tt^tt'^, K~tt~^). The remaining background is some 
combination of non-charm and poorly reconstructed charm events. This is modeled 
with functional forms in mass aexp (6m), energy release aq^^^ + /3q^^^ , and lifetime 
exp(— t/ri) + T] ex.p {—t / T2) ■ Penalty terms are added to ensure backgrounds are 
consistent with known branching ratios. The aS'ected backgrounds are the de- 
cays to K~^K~, 7r+7r~, 7r"'"7r~7r°, ii'^Tr+Tr", Jr~7r+7r°, and K~i~^u which all appear 
in RS and WS. The background due to double misid of the RS signal is fixed based 
on the RS yield and relative efficiency. 

Mini Monte Carlo tests of the fit were used to verify the accuracy of the reported 
fit errors and to check for possible biases in the fit. No significant bias was found 
and the reported fit errors agree with the mini Monte Carlo results. Fit variants 
with different binning, different constraints, and different accoimting of the random 
background were tried. No significant differences were observed. Variations of all the 
selection criteria were also analyzed with no significant differences in the results. For 
the branching ratios reported, systematic errors due to fit variants and cut variants 
were obtained from the r.m.s of the variations and then added in quadrature. 
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Fig. 1. The top plots show projections onto M{D^) and Q{D*) of the fit and data for WS events. 
The bottom-left plot shows the contributions to the WS signal versus t{D°). The bottom-right 
plots shows the 95% CL x' ,y' contour with statistical (inner) and full errors (outer). 



Fitting without mixing we find Rws = (0-4301 o;ogi ± 0.031)%. From the stan- 
dard mixing fit wc find a DCS branching ratio of Rd = (0.382 1 ° ± 0.069)%. 
Figure^shows the WS projections onto M(D°) and Q{D*). The RS data is domi- 
nated by 54452 ± 242 signal events. Due to significant correlations between x' and 
y', the interesting physics result is a contour in the x' , y' plane. Although we fit for 
x'^, we choose to plot x' . The 95% CL contour is defined as the location in the x' , y' 
plane where Alog£ — 2.996 relative to the minimum — log£ in the physical region 
of the x' , y' plane. While the true minimum occurs at x' = —0.06% and y' — 1.0%, 
the minimum with x'^ > occurs at x'^ = and y' — 0.5% with a change in — log>C 
of only 0.006 relative to the true minimum. Systematic checks were performed with 
120 fit and cut variants. The contour variation is consistent with differences in the 
returned value of x' and y' . For the systematic error, we first find the change in 
— log C between the global minimum and the x' , y' location for each of the variants. 
We then find the value greater than 95% of these differences which is 0.482. We find 
the contour at which the change in — log£ is 2.996 -t- 0.482 — 3.478 and call this 
the 95% CL including systematic error. The t{D'~') projection and x',y' contours 
are shown in Fig. ^ Defining 95% CL limits on x' and y' based on the projection 
of the contour onto the respective axis we find x' < 0.83% and 7.2% < y' < 4.1%. 



